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Abstract 

We consider the spatially inhomogeneous and anisotropic reaction- 
diffusion equation u t = m(x)~ 1 div[m(x)a p (x, Vu)] + s~ 2 f(u), involv- 
ing a small parameter e > and a bistable nonlinear term whose 
stable equilibria are and 1. We use a Finsler metric related to the 
anisotropic diffusion term and work in relative geometry. We prove a 
weak comparison principle and perform an analysis of both the genera- 
tion and the motion of interfaces. More precisely, we show that, within 
the time scale of order e 2 | lne|, the unique weak solution u s develops a 
steep transition layer that separates the regions {u e ps 0} and {u £ « 1}. 
Then, on a much slower time scale, the layer starts to propagate. Con- 
sequently, as e — ► 0, the solution u £ converges almost everywhere to 

*AMS Subject Classifications: 35K55, 35K57, 35B25, 35R35, 58B20. 
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in f2^~ and 1 in Q, jy , where and f2j are sub-domains of f2 separated 
by an interface T t , whose motion is driven by its anisotropic mean cur- 
vature. We also prove that the thickness of the transition layer is of 
order e. 

Key Words: nonlinear PDE, anisotropic diffusion, bistable reaction, 
inhomogeneity, singular perturbation, Finsler metric, generation of in- 
terface, motion by anisotropic mean curvature. 



1 Introduction 

Evolution laws for interfaces frequently appear in materials science, differen- 
tial geometry and image processing. In this paper we relate so called diffuse 
and sharp interface models in which interfaces evolve according to an evolu- 
tion law, which involves anisotropic and inhomogeneous driving forces. The 
evolution equations we will consider in particular decrease an inhomoge- 
neous interfacial energy. A diffuse interface model is based on a free energy 
which includes gradient terms and, in this paper, the energy is assumed to 
be of the following Ginzburg-Landau type 

T(u)= [ \a(x,Vu) + \w(u)]m(x)dx , 
Jn ' £ 2 

where Q C R , N > 2, is a bounded domain with smooth boundary, e > 
is a small parameter related to the thickness of a diffuse interfacial layer, 
W is a double well potential with wells of equal depth and m is a positive 
function. We will allow a to be x-dependent and anisotropic, i.e. the value 
of a will depend on the direction of Vu. Taking the gradient flow of T with 
respect to the weighted L 2 -inner-product (u, v) = f n u(x)v(x)m(x)dx leads 
to the following initial boundary value problem for an inhomogeneous and 
anisotropic Allen-Cahn equation 

+ \f(u) infix(0,+oo), 

on dQ x (0, +00), 
in Q, 

where f(u) = — W'(u), v is the Euclidean unit normal vector exterior to d£l 
and a p refers to differentiation with respect to the variable corresponding to 
Vu. We easily derive that solutions to (P £ ) fulfill 

d f 

— f(u) = — {ut) 2 m{x)dx < , 
dt Jn 

i.e. T serves as a Lyapunov function. 



(P £ ) 



I 



Ut 



■ div 



m(x) 
a p (x, Vu) • v = 
u(x,0) = u (x) 



m(x)a p (x, Vu) 
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It can be shown that under appropriate assumptions, see [10], [22], [23] . 
the energies eT converge in the sense of T-convergence to an anisotropic 
functional defined for hyper surf aces, i.e. in the limit e — > the interface is 
sharp. For a smooth hypersurface T the limiting energy becomes 



v / 2a(x, n{x)) m{x)dH N - 1 {x) , 



where n is a suitable Euclidean unit normal vector to V and dTl N ~ 1 refers 
to integration with respect to the (N — l)-dimensional Hausdorff measure. 
Anisotropic energies for hypersurfaces can be analyzed in the context of 
Finsler geometry. If one considers the steepest descent of the anisotropic 
surface energy in relative geometry, where geometric quantities such as cur- 
vature and normal velocity are computed within the context of a Finsler met- 
ric, one obtains an anisotropic and inhomogeneous generalization of mean 
curvature flow. In fact the moving interface Tt evolves according to the law 



(P°) 



m(x) 



y / 2a(x, n) 



Vn 



div 



m{x) 



y / 2a(x,n] 



a p (x, n) 



on r t , 



t=o 



r , 



where V n is the normal velocity of Tt- We will show below that this equation 
can be rewritten in the relative geometry associated with a Finsler metric; 
then it has the form 



(P°) 



on r. 



t=0 



where n^, V n ^ and are, respectively, the anisotropic unit normal in the 
exterior direction, the anisotropic normal velocity of Tt in the direction, 
and the anisotropic mean curvature at each point of Tt- In the isotropic 
homogeneous case one recovers the mean curvature flow V n = —k. We refer 
to a paper by Bellettini and Paolini [5] and Section [2] for details. 

It is the goal of this paper to rigorously prove that Problem (P £ ) con- 
verges to the anisotropic inhomogeneous mean curvature flow (P°), as e — > 0, 
and to give an optimal error estimate between the solutions of (P £ ) and 
those of (P°). We remark that a formal derivation is already contained in 
the paper by Bellettini and Paolini [5]. 

Before going into the details, we note that (P e ) includes the following 
equations as special cases: the spatially inhomogeneous diffusion equation 



ut = dw(A(x)Vu) + \f(u 



(1.1) 
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where A{x) is a positive definite symmetric matrix depending on x; the fully 
anisotropic equation 

ut = div (a p (Vu)) + ^/(u). (1.2) 

The significant difference between (II. ip and (|l,2p is that the anisotropy in 
()1.2p depends on the solution u itself, while it does not in (II. ip . In other 
words, in (|1.2|) the dependence of the energy density on the spatial orien- 
tation of the interface can be chosen much more general when compared 
with (jl.ip where only ellipsoidal energy densities appear. We refer to Gar- 
cke, Nestler and Stoth [15], Barrett, Garcke and Niirnberg [21 [3] for possible 
anisotropic energy densities. Note also that WG allow Clpp (p) to be discontin- 
uous at p = (see Remark [L3] below). 

We suppose in what follows that W(u) is a double-well potential with 
equal well-depth, taking its global minimum value at u = and u = 1. More 
precisely we assume that / = —W is smooth and has exactly three zeros 
< a < 1 such that 

/'(0)<0, /'(o)>0, /'(1)<0, (1.3) 

and that 

/ f(u)du = 0. (1.4) 
J o 

Remark 1.1. Note that we could also consider the case where / is slightly 
unbalanced by order e so that J* f(u)du = 0(e) stands instead of (|1.4j) . 
In this case, the singular limit of (.P £ ) will have an additional driving force 
term in (P°). See Remark 13.21 for details. □ 

The assumptions concerning the anisotropic term are the following. 

(i) a(x,p) is a real valued function, of class Cf^ (for some < ■& < 1) 
on Tlx R N \{0}; 

(ii) a(:r,p) is positive on O x M Ar \{0}; 

(iii) a(x, ■) is strictly convex for all x € fi; 

(iv) a(x,p) is homogeneous of degree two in the p variable, i.e. 

a(x,ap) = a 2 a(x,p) for all (x,p) 6fix R Ar \{0}, all (1.5) 

If p is given by p = (pi, • • • ,£>jv), the vector valued function a p is defined 
by a p (x,p) = (fc(x,p), ■ ■ • , ^(x,p)) , and the matrix valued function a pp 

by a pp (x,p) = ( d f.Q Pi {x,p)\. Moreover, for a vector p = (pi,-- - ,p N ) and 
a matrix ^4 = (ay), we use the notations 

|p|=max|pj| and |A| = max \a,ij\. 

i i,j 
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Remark 1.2. The fact that a is homogeneous of degree two implies that, for 
all (x,p) G U x R^UO}, all a/0, 

o p (x,ap) = aa p (x,p), a pp (x,ap) = a pp (x,p), 
a p (x, ap) ■ p = 2aa(x,p), 
a pp (x,ap)p = a p (x,p). 

By setting a(x,0) = and a p (x, 0) = 0, one can understand that a(x,p) is 
of class C 1 on the whole of x R N . □ 

Remark 1.3. In many important applications in physics, a pp (x,p) is discon- 
tinuous at p = and this makes Problem (P £ ) singularly parabolic. Because 
of lack of uniform parabolicity, our analysis becomes more involved than the 
case (jl.ip or the case of isotropic Allen-Calm equation studied in [I], [llj . 

P2J, HSH23- □ 

We assume that m : $7 — > (0, +oo) is a function of class C 2 such that 
< mi < ?n(x) < ?tt,2 < +oo for any lEll, and that Vm and D 2 m are in 

where D 2 ?n(x) := f ^ ^ (a;)) . 

\dxjdxi ) 

Remark 1.4. Observe that 
1 



m(x) 



div 



m(x)a p (x, Vit) = divo p (x, Vn) + V log m(x) • o p (x, Vw) , (1.6) 



so that our equation contains the generalized Allen-Calm equations discussed 
in Bellettini, Paolini [5], Bellettini, Paolini and Venturini [BJ. □ 

We also assume that the initial data uq G C 2 (r2), and define Co as 

Co := IMI c o(fi) + l|V«o|| C o ( n) + ll^ 2 "o|| c o(n)- ( L7 ) 
Furthermore we define the "initial interface" To by 

To := {x G fi, no(x) = a}, 

and suppose that To is a C 3+,? closed hypersurface without boundary (0 < 
•d < 1), such that, n being the Euclidian unit normal vector exterior to To, 

r CC ft and Vu (x) / if x G r , (1.8) 

uo > a in Qq , «o < a m ' (1-9) 

where Qq denotes the region enclosed by To and the region enclosed 
between d£l and To- 

For T > 0, we set Qt = U x (0,T). We define below a notion of weak 
solutions of Problem {P £ ). For this definition, it is sufficient to only suppose 
that n G H 1 (Q) n L°° (Q) . 
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Definition 1.5. A function u £ G L 2 (0, T; H 1 ^)) n L°°(Q T ) is a weak solu- 
tion of Problem (.P £ ), if 

(i) u% G L 2 (Q T ), 

(ii) a p (^Vti e (3;,t)) G L°°(0, T; L 2 (ft)), 

(iii) it £ (x,0) = uo(^) for almost all x G f2, 

(iv) u e satisfies the integral equality 



for all nonnegative function tp G L 2 (0, T; n L°°(Q T ) and for all 

t€ (0,T). 



One may prove, using monotonicity and compactness arguments as is 
done in [7], [9], that Problem (-P e ) possesses a unique weak solution which 
we denote by u £ . As e — > 0, the qualitative behavior of this solution is the 
following. In the very early stage, the anisotropic diffusion term is negligible 
compared with the reaction term e~ 2 f(u). Hence, rescaling time by r = 
t/e 2 , the equation is well approximated by the ordinary differential equation 
u T = f(u). In view of the bistable nature of /, u £ quickly approaches the 
values or 1, the stable equilibria of the ODE, and an interface is formed 
between the regions {u £ ~ 0} and {u £ ~ 1}. Once such an interface has been 
developed, the anisotropic diffusion term becomes large near the interface, 
and comes to balance with the reaction term so that the interface starts to 
propagate, on a much slower time scale. 

To understand such interfacial behavior, we have to study the singular 
limit of (P £ ) as e — ► 0. Then the limit solution u(x,t) is a step function 
taking the values and 1 on the sides of the moving interface Tt- In the case 
of the usual Allen-Cahn equation, it is well known that Tt evolves according 
to the mean curvature flow V n = — k and we will show in this paper that 
the sharp interface limit of (P £ ) is given by (P°). 

Using the theory of analytic semigroups (see e.g. Lunardi [IB]) it is 
possible to show that the limit Problem (P°) possesses locally in time a 
unique smooth solution. More precisely, there exists a T > such that 
Problem (P°) has a unique solution T = Uo<i<T(^* x M) which satisfies 
r G (7 3 + ,? :( 3 +' ,? )/ 2 . For proofs of the local in time existence of solutions of 
related limit problems, we also refer the reader to [TT] and the discussion at 
the end of Chapter 1 in [16] . 

For each t G (0, T), we define f2j~ as the region enclosed by the hyper- 
surface Tt and Qf as the region lying between <9f2 and Tt- Furthermore we 




(1.10) 
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define a step function u(x,t) by 

I 1 in Qf , , 

^ t ) = { n ■ o- fortG(0,T). (1.11) 

It is convenient to present our main result, Theorem [LBJ in the form of a 
convergence theorem, mixing generation and propagation. It describes the 
profile of the solution after a very short initial period. It asserts that: given 
a virtually arbitrary initial data uq, the solution u £ quickly becomes close 
to or 1, except in a small neighborhood of the initial interface Tq, creating 
a steep transition layer around Tq (generation of interface). The time t e for 
the generation of interface is of order e 2 \ lne|. The theorem then states that 
the solution u e remains close to the step function u on the time interval 
(t e ,T) (motion of interface) . Moreover, as is clear from the estimates in the 
theorem, the thickness of the transition layer is of order e. 

Theorem 1.6 (Generation, motion and thickness of interface). Let n be an 
arbitrary constant satisfying < i] < min(a, 1 — a) and set 

H = /'(a). 

Then there exist positive constants £q an d C such that, for all e G (0, Eq) 
and for almost all (x, t) such that t £ < t <T, where 



we have, 



u £ (x,t) e < 



t e := /_r 1 E 2 \ In e\, (1.12) 

[-77,1+7?] if xeNce^t) 

[-7?, r,] if xefl t -\^ £ (r f ) (1.13) 

[1-7?, 1 + 7?] if xent\M C e{Ttl 



where Af r (Tt) '■= {x S f2, distj,(x, Tf) < r} denotes the r -neighborhood of 
T t ; by dist<f,(x,T t ), we mean the 5$ distance to the set T t , where 5$ is the 
distance associated to a Finsler metric, whose definition is to be given in 
Section [B 



Corollary 1.7 (Convergence). As e — ► 0, the solution u e converges to u 

Jo<i<T(^t 



almost everywhere in Un<t<T(^t = x {*})• 



The organization of this paper is as follows. In Section [21 we recall nota- 
tions and results concerning Finsler metrics that give a natural and efficient 
framework for dealing with anisotropic problems. In Section El we perform 
formal asymptotic expansions in order to derive the equation for the motion 
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of interface, and collect useful estimates on stationary solutions of related 
problems. In Section H] we prove a weak comparison principle for Problem 
(P e ). Such a comparison principle is rather standard, but, in view of the 
fact that a pp (x,p) does not exist at p = 0, we give a complete proof for 
the convenience of the reader. In Section [5l we prove results on the gener- 
ation of interface. For the study of this early time range we construct sub- 
and super-solutions by modifying the solution of the ordinary differential 
equation it t = e~ 2 f(u). In Section El we construct another pair of sub- and 
super-solutions by using the first two terms of the formal asymptotic expan- 
sion given in Section [3l They are used to study the motion of interface in 
the later stage. In Section [TJ by fitting these two pairs of sub- and super- 
solutions together, we prove our main results for (P e ): Theorem 1 1 . 61 and its 
corollary. 

Let us mention some earlier works on anisotropic problems related to 
(P £ ). In [3], Bellettini, Colli Franzone and Paolini study a problem that 
is slightly more general than (P £ ) — by allowing / to be unbalanced in 
the same way as in Remarks 11.11 13.21 of the present paper — and derive a 
very fine error estimate between the formal asymptotic and actual solutions 
of (P e ). We also refer to the articles [EH E], by Elliott and Schatzle, 
on a similar but slightly different problem where the potential W(u) is a 
double obstacle type, namely W(u) = +oo for u £ (0, 1). For the spatially 
homogeneous case a(x,p) = a{p) they prove convergence of the anisotropic 
diffusion problem to an anisotropic curvature flow similar to (P°). Note 
that their second paper [13] considers a kinetic term of the form f3(Vu)ut, 
which makes the meaning of solutions very weak, hence they are treated in 
the framework of viscosity solutions. 

However, in these papers, the authors consider only a very restricted class 
of initial data, namely those having a specific profile with well-developed 
transition layer. More precisely they prove that if the initial data is very close 
to the typical profile that appears in the formal asymptotic expansions of the 
moving interface, then the solution remains close to the formal asymptotic 
for < t < T. In other words the generation of interface from arbitrary 
initial data is not studied there. Summarizing, they have obtained a very 
fine error estimate — of order 0(e 2 ) or higher — between the solutions 
of specific initial data and formal asymptotic, while, in the present paper, 
we consider convergence of solutions of (P £ ) with virtually arbitrary initial 
data to solutions of (P°), with an error estimate of order 0(e). Therefore, 
the two results are both for the convergence of (P £ ) to (P°), but they are 
of different nature. Note that, as far as the thickness of the interface is 
concerned, our 0(e) estimate is optimal (see pQ for details). 

In [8], Benes, Hilhorst and Weidenfeld study both the generation and the 
motion of interface for an anisotropic Allen-Cahn equation which is related 
to ours. Nevertheless, their equation is slightly less general since they do 
not allow x-dependence in a(x,p). Moreover, with their sub- and super- 
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solutions, they cannot achieve the optimal 0(e) estimate of the thickness of 
the interface. 

For numerical simulations for problems (P 6 ) and (P°) we refer to Benes, 
Mikula [9], Garcke, Nestler, Stoth [15], Barrett, Garcke, Niirnberg [21 [3] and 
Paolini [21]. 

2 Finsler metrics and the anisotropic context 

In this section we explain the technique of Bellettini, Paolini [5], Bellettini, 
Paolini and Venturini [BJ to apply Finsler metric to analyze anisotropic non- 
linear problems. The idea is to endow W N with the distance obtained by 
integrating the Finsler metric which makes otherwise lengthy computations 
remarkably simpler. For the convenience of the reader, we first recall basic 
properties of Finsler metrics as stated in [5], [BJ. For more details and proofs, 
see these references. 

2.1 Finsler metrics 

Suppose that (j> : f2 x M. N — > [0, +oo) is a continuous function satisfying the 
properties 

<f)(x,aQ = |a#(x,£) for all (x,£) e!lx R N and all oEl, (2.1) 
Aolfl < <K^,£) < AolCl for all (x,£) € f2 x l w , (2.2) 

for two suitable constants < Ao < Ao < +oo. We say that <p is strictly 
convex if, for any x € CI, the map £ i— > </> 2 (x, £) is strictly convex on M. N . We 
shall indicate by 

B 4> (x) = {^€R N ,<f,(x^)<l} 

the unit sphere of at x G f2. 

The dual function cjp : Cl x — > [0, +oo) of (f> is defined by 

0°(x,r)= sup {£*•£, UB^x)}, (2.3) 

for any (x,£) € x 1^. One can prove that c/> is continuous, convex, 
satisfies properties (|2.1|) and (|2.2|) . and that </> 00 , the dual function of 0°, 
coincides with the convex envelope of <j> with respect to £. 

We say that ^ is a (strictly convex smooth) Finsler metric, and we shall 
write 4> G M(Cl) if, in addition to properties (12. ip and (|2.2I) . ^ and (/>° are 
strictly convex and of class C 2 on SI x W N \ {0}. In particular <ft 00 = <p. 

We denote by 5$ the integrated distance associated to € j\4(Cl), that 
is, for any (x, y) G f2, we set 

<^(x, y) = inf { ^ </>(rf(t)Mt))dt 5 7 G W^ 1 ' 1 ([0, 1]; Cl) , 7 (0) = x, 7 (1) = y}. 

(2.4) 
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In the special case of the Euclidian metric, the function <j> is given by 
(f)(x,p) = (f)(p) = (pi 2 + ••• + pat 2 ) 1 / 2 , so that 5^ reduces to the usual 
distance. 

Given 4> € M(Sl) and sell, let T°(x, •) : R N — > R N be the map defined 

by 

T(x ' u "\o ife = o. ( 5) 

Here (ft® denotes the gradient with respect to p whenever we regard (ft°(x,p) 
as a function of two variables x and p. 

If u : O — > R is a smooth function with non- vanishing gradient, we define 
the anisotropic gradient by 

V^u = T°(x, Vu) = <f>°(x, Vu)^(x, Vu). (2.6) 

If 7] : Q — > is a smooth vector field, we define the m-divergence operator 
by 

div m r] = — — div Imfx)^] = div rt + V log mix) ■ rj , (2-7) 
m(x) 

and then the m-anisotropic Laplacian by 

A^ m tt = div m V^u. (2.8) 

Note that in [S], [B] m is related to while in the present paper m is a 
given function independent of (ft. Nonetheless, in the sequel, we shall use the 
simpler notation A^ := A^ m . 

As in the isotropic case, if Tt is a smooth hypersurface of O at time t, 
and n the outer normal vector to Tt (in the Euclidian sense), we define 
the ^-normal vector to Tt and the 0-mean curvature of Tt by 

n<£ = <f>p(x, n), = div m n$. (2.9) 

Furthermore, if ^ is a smooth function with non-vanishing gradient such 
that T t = {x E Q, ip(x,t) = 0}, and ift is positive outside T t and negative 
inside, then 

n =Sfy n^ = ^(x,VV), (2.10) 

« = div-j^r, = div m 0°(x, V^), (2.11) 

on Tt- We also define the normal velocity of Tt and the ^-normal velocity 
of r* by ^ 

Vn = ~fW = ~<P{.x,vtpy (2 " 12) 

To conclude these preliminaries, we quote a theorem proved in [6]. 
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Theorem 2.1. Let Q be connected, and let 4> G A4(£l). Let 5^ be the 

integrated distance associated to (p. Let C C be a closed set, and let 
dist^^XjC) be the 5^ distance to the set C defined by 

dist rP {x,C)=mf{54x,y) , yeC}. (2.13) 

Then 

4>°(x, Vdistt(x,C)) = 1, (2.14) 
at each point x G £1 \ C where dist^(-, C) is differentiable. 

In the special case of the Euclidian metric, note that (12.141) reduces to 
the property that |Vd| = 1. 

2.2 Application to the anisotropic Allen-Cahn equation 

We set, for all (x,p) G Q x R N , 



</>°(x,p) = ^J2a{x,p). (2.15) 

First, since a(x, •) is homogeneous of degree two, <fp satisfies assumptions 
(pTT|) and {22]) with the constants 

A = [2 min a(x,p)\ 1 / 2 > and A = [2 max a{x,p)\ 1 / 2 > 0. 

a;£fi,|p|=l x£f2,|p|=l 

(2.16) 

By the hypotheses on a(x,p), we see that <fi is strictly convex and of class 
C 2 on fl x R N \ {0}; moreover, by Remark 1 1.2 1 <p° is continuous on the whole 
of n x R N . It follows that 6 is a Finsler metric and the above theory applies. 
We have 

T°(x,p) = { p{ ' W / U ' (2.17) 

if p = 0. 

Let r = Uo<<<t(^< x {*}) be- the unique solution of the limit geometric 
motion Problem (P°) and let d be the signed distance function to T defined 
by 

~ , dist(x,r+) for x G fit, 

d(x,t) = { y ' ; *J (2.18) 

I — dist(x, Ft) for x G Q t , 

where dist(x,I\) is the distance from x to the hypersurface F t in 0. Let d<f, 
be the anisotropic signed distance function to F denned by 

~ . . dist^fa;, I\) for x G fit, 

dJx,t) = { ' ' '_' (2.19) 

I — dist^(z, T t ) for x G 4 , 

where dist<^(x, r 4 ) denotes the 5^, distance to the set F t defined in (|2. 13|) . By 
Theorem 12.11 the following equality holds 

2a{x,Vd < f ) {x,t)) = 1 in a neighborhood of IV (2.20) 
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By setting ip = d and ip = in the second equalities in (|2.10p . (I2.11|) . 
(|2.12p . we obtain two equivalent expressions of the </>-normal vector, the 
0-mean curvature and the ^-normal velocity: 



2a{x,Vd) 



a p (x,Vd) = a p (x,Vd lj) ), 



= div r 



Vn 



1 



2a(x, Vd) 
1 



a p (x, Vd) 



div T 



Op (a;, Vdff,) 



2a(x,Vd) 



(2.21) 

(2.22) 
(2.23) 



The end of this section is devoted to the anisotropic Laplacian 
1 



A^it 



div [m(x)a p (x, Vit) 



m(x) 

divo p (x, Vu) + V log m(x) ■ a p (x, Vit). 



(2.24) 
(2.25) 



In the case of Finsler metrics, it turns out that the term A^u may be less 
regular than Ait. Nevertheless, we show below a boundedness property of 
the anisotropic Laplacian (see [8] for a related property). 

Lemma 2.2. There exists a positive constant Cl such that, for all u E 
C 2,1 (Qt), the following inequality holds: 

[A0u(z,t)| < C L (\Vu(x,t)\ + \D 2 u{x,t)\) for all (x,t) 6 Qt- (2.26) 



Proof. In view of (|2.25p . it is sufficient to deal with the term diva p (x, Vit). 
We can, with no loss of generality, ignore the dependence on time. 

First, assume that x is such that Vu(x) ^ 0. Regarding a(x,p) as a func- 
tion of two variables x and p = (pi, - ■ ■ , p n ), we obtain, by a straightforward 
calculation, that 

,. , _ , d 2 a . . 9 2 a . . <9 2 it . . 

diva p (x, Vu(z)) = 2^ dx . dp . ( x > Vu (x))+2_^ -^r( x > Vtt W)or( 1 )' 



dpidpj 



dxidxj 



(2.27) 



It follows from the homogeneity properties that 



diva p (x, Vu(x))| <|Vit(x)| > max 

yen,|p|=i 



r9 2 



dxjdpj 



(y,p) 



+ \D 2 u(x)\y^ max 



9 2 a 



dpidp. 



■(y,p) 
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where we have used that a is of class C 2 on the compact set Q x = 1}. 
This proves (|2.26p under the assumption Vit(x) 7^ 0. 

Now assume that x is such that Vu(x) = 0. We have to proceed in 
a slightly different way since a pp (x,0) does not make sense. The operator 
a p (x, •) is homogeneous of degree one so that, for any direction £, 

t~ l (a p {x, t() - a p (x, 0)) = a p (x, (). 

We denote by (ei, • • • , e^v) the Euclidian basis of R . It follows from the 
above equality that a p (x, •) admits at the point partial derivatives in any 
direction and 

da p (x, ■) 



dpi 

which, in turn, implies that 



-(0) =aJx, ei ), (2.28) 



±^L( x ,0) = ^L(x,e i ). (2.29) 

Opi Opj Opj 

Note that, since a p (x, •) is homogeneous of degree one, the first term in 
(pT2Tj) vanishes at the point (x,0). It follows from (pT2Tj) and (pT29|) that, in 
the case where Vtt(x) = 0, 

I dl va p (x, Vn(x))| = \^ W (x, ^9^^ 

i,j J J 

da 

max 

. 7/PO.IdI=1 
M 

which proves (|2.26|) in this case as well. □ 



<\D 2 u(x)\Y] max -^-(y,p) 



3 Formal derivation of the interface motion equa- 
tion 

In this section we derive the equation of interface motion corresponding 
to Problem (P £ ) by using a formal asymptotic expansion. The resulting 
interface equation can be regarded as the singular limit of (P £ ) as e — > 0. 
Our argument goes basically along the same lines with the formal derivation 
given by Nakamura, Matano, Hilhorst and Schatzle [21]: the first two terms 
of the asymptotic expansion determine the interface equation. Though our 
analysis in this section is for the most part formal, the results we obtain will 
help the rigorous analysis in later sections. 

Let u e be the solution of (P £ ). Let T = IJo<t<T x W De the solution of 
the limit geometric motion problem and d$ the anisotropic signed distance 
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function to T defined in (I2.19p . We then define 

Q+= |J n+x{t}, Qt= \J ^7><{t}- 

0<t<T 0<t<T 

We also assume that the solution u £ has — one the one hand — the 
outer expansions (away from the interface V), 

u e (x, t) = u(x, t) + euf (x, t) + e 2 uf(x, t) + ■ ■ ■ in Q^, (3.1) 

where u is the step function defined in (II. lip , and — on the other hand — 
the inner expansion (near T) 

u £ {x, t) = U (x, t, + eU^x, t, + e 2 U 2 {x, t, £) + • • • , (3.2) 

near T (the inner expansion), where Uj(x, t, z), j = 0, 1, 2, • • • , are defined 
for x S $7, t > 0, z £ M. The stretched space variable £ := d^(x,t)/e gives 
exactly the right spatial scaling to describe the sharp transition between the 
regions {u £ ~ 0} and {u £ ~ 1}. We normalize Uq in such a way that 

Uo(x, t, 0) = a 

(normalization conditions). To make the inner and outer expansions consis- 
tent, we require that 

U (x,t,+oo) = 1) U k (x,t, +oo) = u£(x,t), 
U (x, t, -oo) = 0, Uk(x, t, -oo) = %(x, t), 

for all k > 1 (matching conditions). 

In what follows we will substitute the inner expansion (|3.2p into the 
parabolic equation in Problem (P £ ) and collect the e~ 2 and e™ 1 terms. For 
this purpose, note that if V = V(x,t,z) and v(x,t) = V(x,t,£) are real 
valued functions then Vv = e~ 1 V z Vd ( f ) + V X V and vt = e~ 1 (d ( j ) )tV z + Vt; 
if v and V are vector valued functions we obtain divv = e~ 1 Vd ( / ) • V z + 
div x V. In the following, we shall use the properties stated in Remark 11.21 
A straightforward computation yields 

u\ = -(d^tUoz + Uot + (fytUu + eU u + --- 

Vu £ = -U 0z Vd^ + V X U + UuVdt + eV x U x + ■■■ 

e 

Op(x, Vu £ ) = -Op(x, UozVdtf, + eV x lI + eUuVd^ + e 2 V x U 1 H ) 

= -a p (x, UozVdfjy) + U 0z Vd^)(V x U + UuVdj) + ■■■ 

Vd i> )(V x U + U lz Vd (f> ) + --- . 
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It follows that 

V log m(x) ■ a p (x, Vu £ ) = -Uq z V log m(x) ■ a p (x, Vd^) + 
and that 

div a p (x, Vu e ) =- Vd^ ■ d z (a p (x, Vu 6 )) + dw x (a p (x, Vu £ )) 



VcL 



Ozz 



a p (x, Vd^) + a pp (x, V d^)^ x U 0z + Ui^Vd^) 

+ - [V x U 0z ■ a p {x, Vd^) + U 0z div a p (x, Vd^)] H 

---^U 0zz 2a(x, Vdtj,) + - 2a(x,Vd,p)U lzz 



+ 2V x U 0z ■ a p (x, Ve^) + U 0z div a p (x, Vd^) 



+ 



where the functions Ui (i = 0, 1), as well as their derivatives, are taken at 
the point (x, £, ^(x, t) / e). Hence, in view of (|2.2U|) . we obtain 

diva p (x, Vu £ ) =^U 0zz + -[U lzz + 2V x U 0z ■ a p (x,Vd^) 

+ U 0z diva p (x, Vd<p)] H . 

We also use the expansion 

f{u £ ) = f(U )+eU 1 f'(U ) + --- . 

Next, we substitute the above expressions in the partial differential equation 
in Problem (P £ ). Collecting the e~ 2 terms yields 

U 0zz + f(U ) = 0. (3.4) 

In view of the normalization and matching conditions, we can now as- 
sert that Uo(x,t, z) = Uo(z), where Uq is the unique solution of the one- 
dimensional stationary problem 

f Uq" + f(U Q ) = 0, 

(3.5) 

{ U (-oo) = 0, U o (0)=a, C/ (+oo) = l. 

This solution represents the first approximation of the profile of a transition 
layer around the interface observed in the stretched coordinates. We recall 
standard estimates on Uq. 

Lemma 3.1. There exist positive constants C and A such that 
< 1 - U (z) < Ce~ xlzl for z > 0, 
< U (z) < Ce~ x \ z \ forz<0. 
In addition to this Uq' > and, for all j = 1, 2, 

\D j U (z)\ < Ce~ xlzl forzeR. 
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Since Uq depends only on the variable z, we have V x Uq = 0. Then, by 
collecting the e _1 terms, we obtain 

U lzz + f'(U )Ux = (d^W - U ' , (3.6) 

which can be seen as a linearized problem for (|3,4|) . The solvability condition 
for the above equation, which is a variant of the Fredholm alternative, plays 
the key role in deriving the equation of interface motion. It is given by 

(d<p) t (x,t) - A^(x,i) U ' 2 (z)dz = 0, 



for all (x,t) € Qt- It follows that (c^)j = A^d^. In virtue of subsection 12.21 
this equation, written in relative geometry, reads as 

V nt< f, = K(f> on F t , (3.7) 

that is the interface motion equation (P°), whereas, in the Euclidian geom- 
etry, the same equation reads as 



m{x) r m(x) 



V n = — div 



y / 2a(x, n) L y / 2a(x^ 



a p (x,n) on T t . (3-8) 



Summarizing, under the assumption that the solution u £ of Problem (P £ 
satisfies 

fl inQ+ 
u — > < _ as e —> 0, 



in Q 



r 



we have formally proved that the boundary Tt between Q t and Qf moves 
according to the law (|3.7|) or ([3. 81) . 

Remark 3.2. To conclude this section, note that (|3.6|) now yields C/i = 0. In 
fact, the second term of the asymptotic expansion vanishes because the two 
stable zeros of the nonlinearity / have "balanced" stability, or more precisely 
because of the assumption J f(u)du = 0. If we perturb the nonlinearity by 
order e, say f(u) < — f{u) — eg(x,t,u), the equation in the free boundary 
problem contains an additional driving force term and U± no longer vanishes. 
More precisely, the equation will read as 

Vn^ = -k^ + cq g(x,t,r)dr onTj, 
J o 

with Co a constant explicitly determined by the nonlinearity /. We refer to 
P] for details. □ 
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4 A comparison principle 



This section is devoted to a comparison principle for weak solutions of Prob- 
lem (P £ ). Such a result is rather standard (see [8]), but, since the problem 
is non-regular where Vit = 0, we prove it here for the self-containedness of 
the paper. 

To begin with, we define a notion of sub- and super-solution of Problem 
(P e ). 

Definition 4.1. A function uf G L 2 (0,T; H 1 (Q)) n L°°{Q T ) is a weak 
super-solution of Problem (P e ), if 

(i) (u+)t G L 2 (Qt), 

(ii) V ( puf(x,t)=a p {x,Vuf(x,t)) G L°°(0, T; L 2 (0)), 
(hi) it £ satisfies the integral inequality 

/ I \(uf) t (p + a p (x,Vuf)-V<p-\f(uf)<p]rn{x)dxdt>0, (4.1) 
Jo Jn 1 el 

for all nonnegative function ip G L 2 (0, T; i? 1 ^)) n L°°(Q T ) and for all 
t€ (0,T). 

We define a weak sub-solution u~ in a similar way, by changing > in 
gU) by <. 

The following remark will be useful when constructing smooth sub- and 
super-solutions in later sections. 

Remark 4.2. If uf G C 2,1 (Qt), it is not difficult to see that uf is a super- 
solution in the sense defined above if and only if 

(i) a p (x, Vuf) ■ v > on d£l x (0, T), 

(ii) Couf > almost everywhere in Qt, 
where the operator Cq is defined by 



C u := u t ]—- div 

m{x) 



m(x)a p (x, Vu) - -^f{u) = u t - A^u - -^f(u). 



In fact, if G C 2,1 (<5t) then, by Lemma 12.21 the function C$uf is well- 
defined in Qt- Also, using Lemma [2T2l we deduce that A^uf G L°°(Qt)- 
The statement is then obtained by integrating (|4.1|) by parts. An analogous 
remark stands for a sub-solution u~ G C 2 ' 1 ^^)- D 

We prove below an inequality which expresses the strong monotonicity 
of the function T°(x,p) = a p (x,p). 
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Lemma 4.3. There exists a constant (3 > such that, for all x € £1, for all 
Pi,P2 G M*, 

(o p (x,p 2 ) - a p (x,pi)) • (p 2 -pi) > /%2 — Pi | 2 - (4.2) 

Proof. First we consider the case that sp\ + (1 — s)p2 ^ for all s E [0, 1]. 
Then, the function s i— > a(x, spi + (1 — s)p2) is of class C 2 on [0, 1] and there 
exist ^3,^4 on the line segment [pi,P2] such that 

a(x,p 2 ) - a(x,pi) = a p (x,pi) ■ (p 2 -pi) + ^(p 2 - pi) • a PP (x,P3)(P2 - pi), 

a(x,pi) - a(x,p 2 ) = a p (x,p 2 ) ■ (pi -p 2 ) + -(pi -P2) • a pp (x,p 4 )(pi -£2)- 

The strict convexity of a(x, •) implies that a pp (x,p) is a positively definite 
symmetric matrix, so that the function (x,p,p) 1— ► a pp (x,p)p ■ p is strictly 
positive and continuous on the compact set Q x S^ -1 x S . Hence there 
exist constants < A 2 < A 2 such that, for all x € fi, all p £ \ {0}, all 

p e M^, 

A 2 |p| 2 <a pp (x,p)p-p< A 2 |p| 2 . (4.3) 

It then follows that 

a(x,p 2 ) - a(x,pi) > a p (x,pi) • (p 2 -pi) + y |p 2 -pi| 2 , (4.4) 

a(x,Pi) - a(x,p 2 ) > a p (x,p 2 ) • (pi -p 2 ) + y |p 2 -pi| 2 . (4.5) 

Adding up inequalities (|4.4p and f|4.5j) yields the desired inequality, with the 
constant (3 = A 2 . 

In the case that sp\ + (1 — s)p 2 = for some s € [0, 1], p\ and p 2 are 
colinear and we may suppose that there exists I € M. such that p 2 = Zpi. 
We can assume Z 7^ 0, I 7^ 1 and pi 7^ 0. By using the properties stated in 
Remark 11.21 we obtain that 

(a p (x,p 2 ) - a p (x,pi)) • (p 2 - pi) = (Z - l) 2 a p (x,pi) - pi 

= 2(Z-l) 2 a(x,pi) 
= 2a(x, (Z - l)pi) 
> Ao 2 |(Z-l)pi| 2 = A 2 |p 2 -pi| 2 , 

where Ao has been defined in (|2.16p . The proof is now completed. □ 

We are now ready to prove the following comparison principle. 

Proposition 4.4 (Comparison principle). Let uf, respectively u~ , be a 
super-solution, respectively a sub-solution, of Problem (P £ ). Assume that 

u 7('ity — U J almost everywhere in fl. 



IS 



Then we have that 

u~ < u £ < uf almost everywhere in Qt- 

Proof. By subtracting inequality (|4.ip for the super-solution uf from in- 
equality for the sub-solution u~ , we obtain that, for all ip G L 2 (0, T; H 1 (Q))n 
L°°(Q T ) such that ip > 0, and for all t G (0,T), 
ft /> 

(«~ — uf)ttp + Vii~) — a p (x, Vuf)) ■ Vip m(x) 

<C f I \u--u+\ip, (4.6) 



JQ, L 



'e "Eiri 
JO Jf7 

where C is a constant depending on e and the L°° norms of /' and m. Next 
we set tp = (u~ - u+)+, which belongs to L 2 (0,T; H 1 ^)) n L°°(Q T ); it 
follows from (|4"T2j) that 

/ (a p (a;, Vu^) — a p (:c, Vm^)) • Vipm(x) 
3 in 

/ Vn~) — a p (x, Vu^)~)) • (Vu~ — Vuf)m(x) 

o J{«7-nt>0} 

>mi/3 [ [ \Vu~ - Vu+| 2 > 0. 

V0 ./{u~-u+>0} 



In view of (14, 6D . we now have that 

* Jo dt Jn K ' Jo J{ u -~u+>o} 

and therefore 

I (<f - «?) + ) 2 (') ^ jf /„ (<«." - «f + ) 2 + /„ (c- - - + > + ) V 

Gronwall's lemma yields 

W - <) + )\t) < e 2Ct/mi J a (W - <) + ) 2 (o). 

Since u~(x,0) < uf(x,0) for almost all x G O, it follows that 

< a.e. in Qt- 

□ 

Lemma 4.5. Let u e be the solution of Problem (P e ) (with initial data uq). 
Then 

-|KIU°°(n) <u £ < max(l, ||ito||z,°°(Q)) a.e. in Q T . 

Proof. By the bistable profile of /, we remark that — ||ito||L°°(Q); respec- 
tively max(l, ||«o||L°°(n))) is a sub-solution, respectively a super-solution, of 
Problem (P £ ). □ 
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5 Generation of the interface 



This section deals with the generation of the interface, namely the rapid 
formation of internal layers that takes place in a neighborhood of Tq = {x £ 
fi, uq(x) = a} within the time span of order e 2 | lne|. In the sequel, rjo will 
stand for the quantity 

rjo := min(a, 1 — a). 
Our main result in this section is the following. 

Theorem 5.1. Let n € (0, t/q) be arbitrary and define [i as the derivative of 
f(u) at the unstable equilibrium u = a, that is 

M = /'(a). (5.1) 
Then there exist positive constants eq and Mq such that, for all e € (0,£o), 

(i) for almost all x £ Q, 

- f] < u £ {x^- l e 2 \\ne\) < 1 + 77, (5.2) 

(ii) for almost all x € ft such that \uq(x) — a\ > Mqe, we have that 

if uo(x)>a + Moe then u £ (x, n~ 1 e 2 \ lne|) > 1 — 77, (5.3) 
if uq(x) < a — Mqe then u £ (x, fi~ 1 e 2 \ lne\) < rj. (5.4) 

We will prove this result by constructing a suitable pair of sub and super- 
solutions. 

5.1 The bistable ordinary differential equation 

The sub- and super-solutions mentioned above will be constructed by mod- 
ifying the solution of the problem without diffusion: 

ut = \f{u), u(x,0) = u (x). 

This solution is written in the form 

u(x,t) = Y (jp2, M^), 

where Y(t, ^) denotes the solution of the ordinary differential equation 

/ n(T,{) =f(Y(T,()) for T>0 

\YM -(. (5 ' 5) 

Here ^ ranges over the interval (— 2Co,2Co), with Cq being the constant 
defined in (jl.7p . We first collect basic properties of Y. 
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Lemma 5.2. We have Y% > 0, for all £ £ (-2C ,2C ) and a// t > 0. 
Furthermore, 

Yc(r = f{Y ^ 



Proof. First, differentiating equation (|5.5p with respect to £, we obtain 

(5.6) 

which can be integrated as follows: 



^(r,£) = exp / 
L Jo 



> 0. 



We then differentiate equation (|5.5p with respect to r and obtain 



y Tr = Y T f(Y), 

Y T (0,0 = f(0, 



which in turn implies 



*r(T,£) = /(£)exp / /'(y( S ,0)^ = /(£)^(r,£). 



(5.7) 



(5i 



(5.9) 



This last equality, in view of (|5.5p , completes the proof of Lemma 15.21 □ 
We define a function A(t,£) by 

/'(Y(r, £))-/'(£) 



/(£) 



(5.10) 



Lemma 5.3. We have, for all £ G (— 2Co, 2Co) and a// r > 0, 

Jo 

Proof. Differentiating the equality of Lemma 15.21 with respect to £ leads to 

Y^ = A(t,OYs, (5.11) 
whereas differentiating (|5.7p with respect to £ yields 



□ 



Yx = Yz / /"(r( s ,£))^( s ,£)d s . 

V o 

These two last results complete the proof of Lemma 15.31 



Next we need some estimates on Y and its derivatives. First, we perform 
some estimates when the initial value £ lies between rj and 1 — 7/. 
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Lemma 5.4. Let rj € (0, 770) be arbitrary. Then there exist positive constants 
C\ = C\(n), Ci = 6*2(77) and C3 = 0-3(7]) such that, for all r > 0, 

(i) if £ E (a, 1 — rj) then, for every r > such that Y(t,£) remains in the 
interval (a, 1 — rj), we have 

Cief^ <Y^r,0 <C 2 e» T , (5.12) 

and 

\A(r,0\<C 3 (e^-l); (5.13) 

(ii) if £ € (77, a) i/ien, /or every r > suc/i £fta£ ^(r, £) remains in the 
interval (jj,a), (I5.12p and (15. 13f) /ioW as u>eZ/, 

where \i is the constant defined in (|5.ip . 

Proof. We take £ € (a, 1 — rj) and suppose that for s G (0,r), y(s,£) 

remains in the interval (a, 1 — rj). Integrating the equality Y T /f(Y) = 1 
from to r yields 



/o f(Y(s,0r" Ji f(q) 
Moreover, the equality of Lemma 15.21 leads to 



ds= I 4r\- ( 5 - 14 ) 



f Y(T ^> f'(q) 



Y ^\m 1 iM^m ]da ( 5i5) 

[ m + m ldq (5 - 15) 



jir + y h(q)dq 



where 

h(q) = (f'(q)-n)/f(q). 

Since 

f W 

the function /i is continuous on [a, 1 — rj]. Hence we can define 

H = H(n) := \\h\\ L oo {ajl _ v) . 

Since \Y(t,£) — £;| takes its values in the interval [0,1 — a — rj] C [0, 1 — a], 
it follows from (|5. 15[) that 

fj,T-H(l-a) < In Y((t,£) < jiT + H{1 - a), 
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which, in turn, proves (|5.12p . Lemma 15, 31 and (|5. 12[) yield 

\A(t,0\ < sup \f"(z)\ [ T C 2 e^ds < C 3 (e" r - 1), 
ze[o,i] Jo 

which completes the proof of (|5.13j) . The case where £ and Y(t, £) are in 
(77, a) is similar and omitted. □ 

Corollary 5.5. Let i] G (0, 770) be arbitrary. Then there exist positive con- 
stants C\ = C\{rj) and C 2 = C 2 (rj) such that, for all r > 0, 

(i) if £, € (a, 1 — 77) then, for every r > such that Y(t, £) remains in the 
interval (a, 1 — 77), we have 

Cxe^H - a) < Y(r,0 - a < C 2 e^ "(£ - a); (5.16) 

(ii) if £ £ (?7, a) i/ien, /or every r > suc/i i/iai l^(r, £) remains in the 
interval (rj,a), we have 

C 2 e» T (£ - a) < y(r,0 - a < C ie ^(£ - a). (5.17) 

Proof. Since 

/(<?)/ (9 - a) -> /'(a) as g -> a, 

it is possible to find B\ = -81(77) > and B 2 = B 2 (rj) > such that, for all 
qe{a,l- r/), 

B 1 (q-a)<f(q)<B 2 (q-a). (5.18) 
We write this inequality for a < Y(t, £) < 1 — 77 to obtain 

B 1 (Y(T,0-a)<f(Y( T ,0)<B 2 (Y(r,0-a). 

We also write this inequality for a < £ < 1 — 77 to obtain 

Bi(£-a)<f(0<B 2 (!i-a). 

Next we use the equality Y^ = f(Y)/f(^) of Lemma [5T2l to deduce that 

^(Y(r,0 - a) < (e - a)y € (r,f) < ^(Y(r,£) - a), 
which, in view of (|5.12|) . implies that 

l^e^ - a) < F(r,e) - a < ^C 2 e^(£ - a). 
t> 2 B\ 

This proves (|5.16p . The proof of (|5.17p is similar and omitted. □ 

Next we present estimates in the case where the initial value £ is smaller 
than 77 or larger than 1 — rj. 
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Lemma 5.6. Let 77 G (0,170) an d M > be arbitrary. Then there exists a 
positive constant C4 = 6*4(77, M) such that 

(i) if £ G [1 — 77, 1+M], then, for allr > 0, Y(t,£) remains in the interval 
[1 - 77, 1 + M] and 



\A(t,0\ <C±t for t >0; 



(5.19) 



(ii) if £ G [— M, 77], then, for all r > 0, Y(t, £) remains in the interval 
[— M, 77] and (15.190 ZioWs as we//. 

Proof. Since the two statements can be treated in the same way, we will 
only prove the former. The fact that Y(r, £), the solution of the ordinary 
differential equation (|5.5p . remains in the interval [1 — 77, 1 + M] directly 
follows from the bistable properties of /, or, more precisely, from the sign 
conditions /(l - 77) > 0, /(l + M) < 0. 

To prove (|5T9j) . suppose first that £ G [1, 1 + M]. In view of CL3|), /' 
is strictly negative in an interval of the form [1,1 + c] and / is negative in 
[1, 00). We denote by — m < the maximum of / on [1 + c, 1 + M]. Then, as 
long as Y(t, £) remains in the interval [1 + c, 1 + M], the ordinary differential 
equation (15. 5p implies 

Y T < —m. 

By integration, this means that, for any £ G [1, 1 + M], we have 

M-c 



Y(t,£) G [1,1 + c] for t >t :-- 



m 



In view of this, and considering that f'(Y) < for Y G [1,1 + c], we see 
from the expression (|5.7p that 



F € (r,0 =exp / /'(Y(s,0)ds 



exp 



/'(y( S ,0)d« 



< 



exp 



f'(Y(s,0)ds 







< exp 



sup =: C 4 = C 4 (M), 

ze[-A/,l+A/] - 1 



for all t > r. It is clear from the same expression (|5.7p that 1^ < C4 holds 
also for < r < r. We can then use Lemma 15.31 to deduce that 

\A(t,0\ <cJ T \f"(Y(s,0)\ds 



< CA sup ze[ _ M>1+M] \f"(z)\\T =: C 4 r 
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The case £ € [1 — 77, 1] can be treated in the same way. This completes the 
proof of the lemma. □ 



Now we choose the constant M in the above lemma sufficiently large 
so that [-2C ,2C ] C [-M, 1 + M], and fix M hereafter. Then C 4 only 
depends on 77. Using the fact that r = 0(e^ T — 1) for r > 0, one can easily 
deduce from f|5. 13j) and (|5.19|) the following general estimate. 

Lemma 5.7. Let 77 € (0, 770) be arbitrary and let Co be the constant defined 
in (|1.7|) . Then there exists a positive constant C5 = 6*5(77) such that, for all 
£ € (-2C , 2C ) and a// r > 0, 

L4(T,0|<C 5 (e" T -l). 



5.2 Construction of sub- and super-solutions 



We are now ready to construct sub- and super-solutions in order to study the 
generation of the interface. By using some cut-off initial data, see subsection 
3.2 in PP, we can modify slightly uq near the boundary dSl and make, without 
loss of generality, the additional assumption 



a p (x, Vuq(x)) ■ v = 
Our sub- and super-solutions are defined by 



on 



wf(x,t) = Y[-=, u (x) ±e 2 C 6 (e 



fit/e 2 



r 



(5.20) 



(5.21) 



Lemma 5.8. There exist positive constants Eq and Cq such that, for all 
e G (0,£o), (w~,wf) is a pair of sub- and super- solutions of Problem (P £ ), 
in the domain Q x (0, [i~ 1 £ 2 \ lne|). 

Proof. Following Remark 14.21 we define the operator Cq by 



C u := u t 



1 



m(x) 



div 



m(x)a p (x, Vn) 



1 



/(«) 



(5.22) 



and prove that Cowf > 0. We compute 

£ z 

Vio+ = Vu (x)Y i . 

Using (|5.2U|) and the fact that a p (x, ■) is homogeneous of degree one, we see 
that wf 1 satisfy the anisotropic Neumann boundary condition a p (x, J \/w^ z ) ■ 
v = on dQ x (0, +00). In view of the ordinary differential equation (|5.5I) . 
we obtain 



C w£ = fiC 6 e^ £ \ 



1 



mix 



■ div 



m(x)a p (x, Vwf) 
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By the estimate of the anisotropic Laplacian (|2.26|) . it follows that 

C wt > ^C^ t l e2 Y i - C L (\Vwf(x,t)\ + \D 2 wf(x,t)\) , (5.23) 

where we recall that \D 2 wf(x, t) \ = maxjj \didjwf(x, t)\. A straightforward 
calculation yields 

didjwf(x,t) = (didjU )Y£ + {diU Q djU Q )Y^. 

Recalling that 1^ > 0, we now combine the expression of Vwj", the above 
expression and inequality (|5.23p to obtain 

Cowt/Y^ > ^C 6 e^ £2 - C L C -C - C 2 ^ , (5.24) 

where Co is the constant defined in (II. 7p . We note that, in the range 
(0, /i _1 e 2 | hie|), we have 

< e 2 C 6 (e^ £2 - 1) < e 2 ^- 1 - 1) < C , 

if £q is small enough. Hence 

£ := u (x) ±e 2 eC 6 (e^ £2 - 1) € (-2C ,2C ), 

so that, by the results of the previous subsection, Y remains in (— 2Cq, 2Cq). 
In view of (|5.11|) . Y^/Y^ is equal to A so that, combining the estimate of A 
in Lemma 15.71 and (|5.24p . we obtain 

C-owt/Y^ > ( M C 6 - C 2 C 5 )e^ £2 - C L C - C . 
Now, choosing 

C 6 > -max(Co 2 C 5 ,Cb(C L + l)) 

proves C^wf /Yjc > 0. Since Y% > 0, it follows that Cowf > 0. Hence, 
by Remark 14.21 w t ls a super-solution of Problem (P £ ). Similarly w~ is a 
sub-solution. Lemma loTHI is proved. □ 

To conclude this subsection, we remark that w^(x, 0) = y(—^, uq(x 
uo(x). Consequently, by the comparison principle, 

w~(x,t) < u £ (x,t) < w£(x,t) , (5.25) 

for almost all (x,t) G SI x (0, fi~ 1 e 2 \\ne\). 
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5.3 Proof of Theorem I5TT1 

In order to prove Theorem 15. II we first present a key estimate on the function 
Y after a time interval of order r ~ | lne|. 

Lemma 5.9. Let rj G (0, r?o) be arbitrary; there exist positive constants So 
and Cj such that, for all e G (0, £$), 

(i) for all £ G (-2C , 2C ), /or a// r > me|, 

-7?<y(r,0 < l + r?, (5.26) 

(ii) for all £ G (— 2Co,2Co) smc/i i/iai |£ — aj > CVe, /or a// r > /j" 1 ) lne|, 

if ^>a + C 7 e then Y(t, f ) > 1 - rj, (5.27) 
if £<a-C 7 £ then F(r,£) < rj. (5.28) 

Proof. We first prove (|5.27[) . For £ > a + Cj£, as long as Y(t, £) has not 
reached 1 — r], we can use (I5.16D to deduce that 

Y(t, £) > a + Cie" T ($ - a) > a + CiCVe^e > 1 - r? , 

provided that r satisfies r > /j" 1 In !+ r "! ? . Choosing 



c 7 



ClC 7 £ 

max(a, 1 — a) — 77 



completes the proof of (|5.27p . Using (I5.17|) . one easily proves (15.281) . 

Next we prove ()5.26j) . First, by the bistable assumptions on /, if we 
leave from a £ G [—7/, 1 + 7/] then Y(t, £) will remain in [—rj, 1 + 77]. Now 
suppose that 1 + 7] < £ < 2C . We check below that l^" 1 ! lne|, < 1 + 77. 
First, in view of (II. 3p . we can find p > such that 



if Ku<2C then /(tt) < p(l - u), 
if - 2C < u < then /(u) > -pu. 

We then use the ordinary differential equation (|5.5[) to obtain, as long as 
1 + rj < Y < 2C , the inequality Y T < p(l - Y). It follows that 

< -P- 



F- 1 

Integrating this inequality from to r leads to 

Y(t, £) < 1 + (£ - l)e" pr < 1 + (2C - l)e" pr . 

One easily checks that, for e G (0,£o) ? with £0 = £o( r ?) small enough, we 
have Y(t, £) < 1 + 77, for all r > /j _1 |lne|, which completes the proof of 
(I5T26]) . □ 
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We are now ready to prove Theorem 15.11 By setting t = fi l e 2 \ lne| in 
(I5.25p . we obtain, for almost all x € 0, 

Y^llneluoix) - (C 6 e - C e e 2 )) 

< u E (x,fi~ 1 £ 2 \lne\) < y(// _1 |lne|,ito(x) + C 6 e-C 6 e 2 ). (5.30) 

Furthermore, by the definition of Co in (jl.7p . we have, for eq small enough, 

-2C < uq(x) ± (C 6 e - C 6 e 2 ) < 2C for xEfl. 

Thus the assertion (|5.2|) of Theorem 15.11 is a direct consequence of (|5.26p 
and (IQnjl . 

Next we prove (|5.3p . We choose Mo large enough so that Mqe — CqE + 
Cq£ 2 > C^e. Then, for any x G Q such that tto(x) > a + Mo£, we have 

uq(x) - (C 6 e - C 6 e 2 ) > a + M e - C 6 e + C 6 e 2 >a + C 7 e. 

Combining this, (I5.30|) and (|5.27p . we see that 



for almost all x G that satisfies uq{x) > a + Mq£. This proves (|5.3p . The 
inequality ()5.4p can be shown the same way. This completes the proof of 



6 Motion of the interface 

We have seen in Section [5] that, after a very short time, the solution u e 
develops a clear transition layer. In the present section, we show that it 
persists and that its law of motion is well approximated by the interface 
equation (P°). 

More precisely, take the first term of the formal asymptotic expansion 
(|3.2p as a formal expansion of the solution: 



The right-hand side of (16. ip is a function having a well-developed transition 
layer, and its interface lies exactly on Tf. We show that this function is a 
very good approximation of the solution; therefore the following holds: 

If u £ becomes rather close to u £ at some time moment, then it 
stays close to u e for the rest of time. 



u 6 (x,fi 1 e 2 \ lne|) > 1 — r\ 



Theorem 15.11 



□ 




(6.1) 
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To that purpose, we will construct a pair of sub- and super-solutions u~ 
and uf of Problem (P £ ) by slightly modifying u £ . It then follows that, if 
the solution u £ satisfies 

uj(x,t ) < u £ (x,t ) < uf(x,t ) , 

for some to > and for almost all x £ fi, then 

uj(x, t) < u e (x, t) < uj(x, t) , 

for almost (x, t) € Qt that satisfies to < t < T. As a result, since both 
uf ,u~ stay close to u £ , the solution u £ also stays close to u £ for to < t < T. 

6.1 Construction of sub and super-solutions 

To begin with we present a mathematical tool which is essential for the 
construction of sub and super-solutions. 

A modified anisotropic signed distance function. Rather than work- 
ing with the anisotropic signed distance function d^, defined in (I2.19p . we 
define a "cut-off anisotropic signed distance function" as follows. Choose 
do > small enough so that g?^(-, •) is smooth in the tubular neighborhood 

of r 

{(x,t) £ Q T , \d^{x,t)\ < 3d }, 

and that 

dist (r t , dfl) > 3d for all t G (0, T). (6.2) 
Next let £(s) be a smooth increasing function on R such that 



COO 



s if \s\ < do 

-2d if s < -2d 
2d if s > 2d . 



We define the cut-off anisotropic signed distance function d^ by 

d < f > (x,t) = ((d cf> {x,t)). (6.3) 

Note that, in view of (pT20j) . 

2a(x, Vd^(x, t)) = 1 in a neighborhood of I\, (6-4) 

more precisely in the region {(x,t) E Qt, |d^(x,t)| < do}- Moreover, in 
view of (|6.2|) . we have 

2a(x, Vd<f>(x, t)) = far away from T t , (6-5) 
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more precisely in the region {(x,t) £ Qt, \d^{x^t)\ > 2do}. Furthermore, 
since the moving interface T satisfies Problem (P°), an alternative equation 
for r is given by 

{d(j,)t = —]—r div [m(x)a p (x, Vd^)] on T t . (6.6) 



Construction. We look for a pair of sub- and super-solutions uf for (P £ ) 
of the form 

u e {x,t) = U ^ J±g(t), (6-7) 

where Uq is the solution of (|3.4p . and where 

p(t) = - e -/»/e a + e L * + K, 
9 (t) = tr (p e -Pt/*? +e 2 Le Lt ). 

Note that g = ae 2 pt- It is clear from the definition of uf- that 



[ 1 for all (x, t) G 
limit? fx, t) = <^ T (6i 

1 for all (z,t) G Q^. 



The main result of this section is the following. 

Lemma 6.1. There exist positive constants (3, a with the following proper- 
ties. For any K > 1, we can find positive constants so and L such that, for 
any e £ (0, Eq), the functions u~ and uf satisfy the anisotropic Neumann 
boundary condition and 

CqU~ < < Cou+, 
in the range O x (0, T), where the operator C$ has been defined in (15.221) . 



6.2 Proof of Lemma 16.11 

We show below that 

C uf := (u+)t TT div [m(x)a p (x, Vu+ )] - \f(uf) > 0, 

m(x) e z 

the proof of inequality Lqu~ < follows by similar arguments. 
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6.2.1 Computation of Couf 

In the sequel, the function Uq and its derivatives are taken at the point 
(d ( f ) (x,t) + ep(t))/e. Straightforward computations yield 

(u+) t = (q(d4,)t+Pt)Uo' + q u 

diva p (x, Vuf) = ^U "Vd^ ■ a p {x,Vd ( f > ) + -Uq div a p (x, Vd^) 

= -^Uo"2a(x, Vd$) + -Uq' div ctp(x, Vd^), 

where we have used properties stated in Remark 11.21 Note that, d^ being 
constant in a neighborhood of d£l, we have that Vu+ = on dtt x (0, T) and 
uf satisfies the anisotropic Neumann boundary condition a p (x, Vu+) • v = 
on dQ, x (0, T). At last, we use the expansion 

f( u t) = f(U )+ q f'(U ) + l q 2 f"(9), 

for some function 6(x,t) satisfying Uq < 9 < uf . 

Combining the above expressions with (|1.6p and (13. 5D , we obtain CquJ = 
E± + £/2 + E% , where 

#1 = -^q{f{Uo) + \qf"{9)) + U 'pt + g t , 

U " / \ 
E 2 = -p-(l-2o(x > Vd^)J > 

^3 = "^-((^)t - ^y div [m{x)a p {x, 

In order to estimate the above terms, we first present some useful in- 
equalities. As /'(0) and /'(l) are strictly negative, we can find strictly 
positive constants b and m such that 

if Uq(z) € [0, b] U [1 - 6, 1] then f'(U {z)) < -m. (6.9) 

On the other hand, since the region {(x, z) £ fl x M| Uq(z) £ [b, 1 — b] } is 
compact and since C/q' > 011 there exists a constant ai > such that 

if Uq{z) € [b, 1 - 6] then f7 '(» > ai. (6.10) 

We now choose M > such that |J7o| < M — 1. We then define 

F= sup + + (6.11) 

\z\<M 
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P = J , (6-12) 

and choose a that satisfies 

< a < min(cro, a\, a 2) , (6.13) 

where 

ai 1 4/3 

ffo: -i^' ai: -^+i' a2 '-F(pTTy 

Hence, combining (16. 9p and (16.101) . we obtain, using that a < oq, 

U '{z) - af'(U {z)) > Aa(3 for z e M. (6.14) 

Now let -ftT > 1 be arbitrary. In what follows we will show that Couf > 
provided that the constants £0 and L are appropriately chosen. From now 
on, we suppose that the following inequality is satisfied: 

elLe LT < 1 . (6.15) 

Then, given any e G (0,£q), since a < a±, we have < q(t) < 1, hence 

- M <uf(x,t) < M . (6.16) 

6.2.2 An estimate for E\ 

A direct computation gives 



e 

where 



E 1 = ^ e-^ 2 (I - aP) + Le Lt (I + e 2 aL), 



I = U ' - vf'(U ) - y f"(e)(Pe~W* 2 + e 2 Le Lt ). 



In virtue of (|6.14p and (I6.16|) . we obtain 



2 



I>4a/3- -j F (P + £ Le 
Then, in view of (|6.15p . using that a < 02, we have I > 2a (3, which implies 



El > + 2af3Le Lt =: ^e'^ + Le Lt . 



32 



6.2.3 An estimate for E2 

First, in the points where where < do, by (16, 4p . we have E2 = 0. Next 
we consider the points where \d^\ > do. We deduce from the definition of Ao 
in (gJSD that 

< 2a(x,Vd (f> (x,t)) < (A°) 2 |Vd (x,t)| 2 

<(A ) 2 ||V<y2,:=.D<cx). 

Applying Lemma 13.11 yields 

I-E2I < -^(l + m e - A l^ +£p l/ £ < ^(l + D)e-^ d °/ e -^l 

We remark that < K — 1 < p < e LT + K, and suppose from now that the 
following assumption holds: 

Then — - |p| > so that, defining C" := C(l + D), 



e 2 ~ (eAd ) 2 ' 

6.2.4 An estimate for £"3 

We set 

G(x, t) = {ds)t{x, t) ]—- div [m(x)a p (x, VdJx, t))] . 

m{x) 

We recall that d^ G C 3+1? ' ( 3+1? )/ 2 in a neighborhood V of T, say 

V = {(x,t) G Q T , |^(M)I < ^o}- 
Combining the fact that 

2a(x,S7d <t> (x,t)) = 1 in V, 
with the definition of A in (|2.16p . we see that 

|V<g > ^ in V. (6.18) 



We also recall that (x,p) 1— > a(x,p) is of class on fi x 1^ \ {0}. Since 

|Vd^| is bounded away from zero, it follows that x 1— > Q{x,t) is Lipschitz 
continuous on V. By equation (|6.6p . we have that 



Q(x,t) =0 on T t = {x G VL,dA){x,t) = 0}, 



33 



and it follows from the mean value theorem applied separately on both sides 
of Ft that there exists a constant N\ such that 

\G(x,t)\ < JVi|cfy(sG,t)| for all (x,t) € V. (6.19) 

Next, using Lemma 12.21 we remark that Q is bounded on x [0,T]\V so 
that there exists a constant N 2 such that 

_ sup \G(x,t)\ < N 2 . (6.20) 

C2x[0,T]\V 

By the inequalities (|6.19|) and (|6.2U|) . we deduce that 
\Q(x,t)\ < N\di(x,t)\ mQ T , 
with N := max(iVi, N 2 /do). Applying Lemma l3~Tl we deduce that 

\E 3 \ < NC W e -A|<Ve+p| 

e 

< NCmax ym \y\e- x \ y+p \ 

< NCmax(\p\ ,-). 

A 

Thus, recalling that \p\ < e Lt + K, we obtain 

\E 3 \<C 3 (e Lt + K)+C 3 ', 
where C 3 := NC and C 3 ' := NC/X. 

6.2.5 Completion of the proof 

Collecting the above estimates of E\, E 2 and E 3 yields 

Cou+ > ^e-^ 2 + (LCV - C 3 )e Lt - C A , 
where C4 := C 2 + KC 3 + C 3 . Now, we set 

T 4e ' 

which, for £0 small enough, validates assumptions (16.150 and ()6.17p . If 
£0 is chosen sufficiently small (i.e. L sufficiently large), we obtain, for all 
£ G (0,£ ), 

£ uf > (LC\' - C 3 )e Lt -C A > -LC ± ' - C 4 > 0. 
The proof of Lemma 16. II is now completed. □ 
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7 Proof of Theorem 11.61 and Corollary 11.7 



(7.2) 



Let r\ G (0, rjo) be arbitrary. Choose (5 and a that satisfy (|6.12p . f)6. 13j) and 

<tf < |. (7.1) 

By the generation of interface Theorem 15. 1\ there exist positive constants 
£o and Mq such that t)5.2[) , (15.30 and (15.4)) hold with the constant r\ replaced 
by cr/3/2. Since, by the hypothesis (|1.8|) and the equality (|2.21|) . Vuo(x) • 
n<f>(x) / everywhere on the initial interface Tq = {x G 0,, uq(x) = a} and 
since Tq is a compact hypersurface, we can find a positive constant M\ such 
that 

if d^(x,0)> M\£ then n (x) > a + M e, 
if d^,(a;, 0) < — M\E then uo( ;c ) < « — Mo£. 
Now we define functions H + (x), H~{x) by 

ff+n = J 1 + (T/9/2 if g^(z,0) > -Mie 
W 1 0^/2 if d^(x,0) < -Mie, 

f 1 - cr/3/2 if ^(x,0) > Mie 
W \ -cj/3/2 if d^(x,0) < Mie. 

Then from the above observation we see that 

H~{x) < u £ (a;,// _1 £ 2 |lne|) < H + (x), (7.3) 

for almost all x £ fi. 

Next we fix a sufficiently large constant K > 1 such that 

[/ (-Mi + #)>l-^ and ^(Mi-tf) < ^. (7.4) 

For this -fiT, we choose eq an d L as in Lemma 16.11 We claim that 

u~(x,0) < H~(x), H + {x) < u+(x,0), (7.5) 

for all x £ £1. We only prove the former inequality, as the proof of the latter 
is virtually the same. Then it amounts to showing that 

u~(x,0) = U (^^ -K) -<j{(3 + e 2 L) < H~{x). (7.6) 

In the range where d<f,(x,0) < M±e, the second inequality in (|7.4j) and the 
fact that Uq is an increasing function imply 

t/ ( ^ (X ' 0) — -fiT) — a(P + e 2 L) < U (M 1 — K) — aft — ae 2 L 

< H~(x). 
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On the other hand, in the range where d^x, 0) > M±£, we have 

< H-(x). 

This proves (|7.6|) . so that (|7.5|) is established. 
Combining (|7.3f) and (|7.5j) . we obtain 

ti~(rr,0) < u £ (x, /U _1 e 2 | lne|) <ii+(x,0), 

for almost all Since, by Lemma 16. 11 u~ and u+ are sub- and super- 

solutions of Problem (P £ ), the comparison principle yields 

u~{x,t) < u £ (x,t + t £ ) < u+(x,t), (7.7) 

for almost all (x,t) € Qt that satisfies < t < T — t £ , where we recall 
that t £ = /i _1 e 2 |lne|. Note that, in view of (|6.8|) . this is sufficient to prove 
Corollary 11.71 Now let C be a positive constant such that 

U (C-e LT -K) > 1-| and U (-C + e LT + K) < ^. (7.8) 

One then easily checks, using (17.70 and (|7.1|) . that, for eo small enough 
and for almost all {x, t) € we have 

if d,j,(x,t)> Ce then u £ {x, t + i e ) > 1 - rj, 
if d ( f ) (x,t) < —Ce then n e (x,t + t £ ) < r], 

and 

-u e (x,t + t e ) € [—77, l+r7], 
which completes the proof of Theorem 11.61 □ 
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